Unitarity bound for gluon shadowing 
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Although at small Bjorken x gluons originated from different nucleons in a nucleus overlap in the 
longitudinal direction, most of them are still well separated in the transverse plane, therefore cannot 
fuse. For this reason the gluon density in nuclei cannot drop at small x below a certain bottom 
bound, which we evaluated in a model independent manner assuming the maximal strength of gluon 
fusion. We also calculated gluon shadowing in the saturated regime using on the Balitsky-Kovchegov 
equation, and found the nuclear ratio to be well above the unitarity bound. The recently updated 
analysis of parton distributions in nuclei [l|] including RHIC data on high-pr hadron production at 
forward rapidities, led to astonishingly strong gluon shadowing, which is far beyond the unitarity 
bound. This indicates a misconception in the interpretation of the nuclear suppression observed at 
HRIC. 

PACS numbers: 24.85.+P, ll.80.La, 13.60.Hb, 13.85.Rm 
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I. INTRODUCTION 

Gluon shadowing, or suppression of the gluon density 
at small Bjorken x in nuclei, keeps challenging the high- 
energy-physics community. Lacking direct information 
from data, attempts have been made to extract gluon 
density from data on deep-inelastic scattering (DIS) on 
nuclei. Parton distribution functions (PDF) for different 
species are related to each other via the DGLAP evolu- 
tion, and the dependence of the quark distribution in 
nuclei should contain information about the gluon den- 
sity. We think that the best analysis of this kind done 
so far was performed in the next-to- leading order (NLO) 
by De Florian and Sassot [1] , who found good sensitivity 
of data to the gluon distribution and concluded with a 
rather weak effect of gluon shadowing (as was predicted 
in Bi). 

However, a word of caution is in order. This kind of 
analysis cannot be considered as sefconsistent from the 
theoretical point of view. Indeed, an NLO analysis has 
to include an equation for gluon shodowing with ac- 
curacy. However, such an equation has not been found 
so far, therefore NLO analyses 0, Q are not reliable. 

A self-consistent theoretical scheme has been formu- 
lated so far only for a leading order (LO) analysis. Un- 
fortunately, previous attempts to obtain the gluon distri- 
bution in nuclei employing a LO procedure f^, ^ failed 
since the data from the New Muon Collaboration (NMC) 
experiment involved in the analyses spans only down to 
X « 10~^ and do not have sufficient accuracy to be sensi- 
tive to the gluon distribution (the frequently used gluon 
shadowing from [6| is just an educated guess). 

In order to reach the smallest values of x in nuclei one 
could rely on data for hard processes in hadron-nucleus 
collisions at forward rapidities. Indeed, the BRAHMS ex- 
periment Q at RHIC performed measurements of high- 



Pt pion production in c? — Au collisions at ^/s = 200 GeV 
and pseudo-rapidities r/ = 2.2 and 77 = 3.2, and the STAR 
experiment reached 77 = 4, in the deuteron fragmen- 
tation region. Although these data allow to access a new 
range of very small x in the nucleus, another word of cau- 
tion is in order: QCD factorization breaks down close to 
the kinematic limit 10]. This happens because the phase 
space for partons in the target does not rise, but shrinks 
(both for proton or nuclear targets), and because of ini- 
tial/final state interactions in the nucleus. Indeed, inclu- 
sion of these data in the new global analysis [l[ resulted in 
an astonishingly small amount of gluons in nuclei. Fig. [T] 
shows the result for the gold to proton ratio of gluon 
densities R^{x,Q^) as function of x at — 1.69 GeV^. 
According to [1] gluons in the nucleus are suppressed by 
orders of magnitude at small x. This implies that the 
effective number of nucleons, 



A, 



i,ff=ARf{x,Q-'), (1) 

is extremely small. At cc = 10~^ the new nuclear PDF 3 
predicts A^ff = 1, i.e. the whole gold nucleus is repre- 
sented by a single nucleon (!). Certainly, this cannot be 
true, unless all nuclear gluons are strongly correlated in 
impact parameter plane throughout the whole nucleus. 
Such a strong gluon shadowing also contradicts the bot- 
tom unitarity bound derived below, which is plotted in 
Fig. [1] by the thin curve. The nuclear ratio resulted from 
the EPS08 analysis is fully in the shaded area forbidden 
by the bound. 

This controversy indicates a misconception in the in- 
terpretation of the suppression effect observed by the 
BRAHMS experiment. This also explains the failure of 
the analysis m to describe data ^Sj at < 4 GeV^ and 
similar data at larger rapidity 77 = 4. 

Actually, the effect of enhanced nuclear suppression 
at forward rapidities observed by the BRAHMS exper- 
iment is not unique. It is known to happen in every 
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FIG. 1: Gold to proton ratio of gluon densities at Q — 
1.69 GeV^ as function of Bjorken x. The thin curve shows 
the bottom unitarity bound for the gluon ratio Rf{x), and 
the shaded area covers the values forbidden by this bound. 
The thick solid curve is calculated with the EPS08 code 



reaction measured so far, even at low energies, where no 
coherent effects are possible (see examples in Al- 
though at forward rapidities one accesses the smallest x 
in the nuclear target, one simultaneously gets into the 
region of large x of the beam hadron, where energy con- 
servation becomes an issue. The projectile hadron, or its 
remnants, propagating through the nucleus dissipate en- 
ergy for multiple interactions in the nucleus. This makes 
it more difficult to give the main fraction of the initial 
energy to one particle detected with large Feynman xp 
or/and transverse xt- Consequently, multiple interac- 
tions enhance nuclear suppression at forward rapidities. 
Calculations performed in 10] based on the AGK cut- 
ting rules well describe the BRAHMS data [8| and cor- 
rectly predict [ll| suppression at larger rapidity as 
well as in other reactions measured at large xp- This 
nuclear suppression scales (approximately) in Feynman 
xp, rather than in X2 as should have been expected for 
nuclear shadowing. Notice that the effective energy loss 
is a nonperturbative effect, which can be evaluated only 
within phenomenological models. 

The effect of energy dissipation caused by multiple col- 
lisions should not be mixed up with gluon shadowing 
which is also related to multiple interactions, but is re- 
sponsible only for a part of the observed nuclear effects. 
Attributing the whole suppression to gluon coherence ef- 
fects (shadowing or color glass) leads to grossly exagger- 
ated gluon shadowing [l| which, as we will show below, 
breaks the unitarity bound. Therefore, the derivation of 
the minimal value for may give a useful guide for the 
separation of the two effects, gluon shadowing and en- 
ergy loss, in the analysis of experimental data. In order 



to discriminate the two mechanisms, one can also move 
to lower energies where coherence phenomena disappear, 
while the effects of energy conservation, which scale in 
xp, remain essentially unchanged. 



II. GLUON SATURATION 

Shadowing, both the term and the phenomenon, which 
came from optics, got a new interpretation within the 
parton model [l2| . Although the bound nucleons are well 
separated in the nuclear rest frame, they start "talking" 
to each other being boosted to the infinite momentum 
frame. At first glance both the inter-nucleon spacing and 
the nucleons themselves are subject to Lorentz contrac- 
tion with the same 7-factor, 7^? = p^/mjv. So they 
should remain separated. However, a boosted nucleon 
develops quantum fluctuations, sea partons, which are 
an analogue to Weizsacker- Williams photons. A parton 
carrying a fraction x of the proton light-cone momentum 
contracts less, since has a smaller 7-factor, 7(2;) — xjn. 
Thus, the parton clouds originated from different nucle- 
ons should overlap at small a; < < 1 and may fuse leading 
to a reduction of the parton density. The correspond- 
ing nonlinear fusion term was introduced in the Gribov- 
Levin-Ryskin [l| and Mueller-Qiu [11] (GLR-MQ) equa- 
tion. 
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Here 



G^(x,Q^6) 



Ga{x,Q^), 



(3) 



is the impact parameter density of the nuclear gluon dis- 
tribution function. In absence of the nonlinear term in 
(PI) the unshadowed nuclear gluon density can be related 
to the nucleon one as, 



G''A{x,Q^b)=TA(b) Gn{x,Q^ 



(4) 



where TA{b) = dz pA{b, z) is the nuclear profile func- 
tion, the integral of the nuclear density along the hadron 
trajectory. Thus, the gluon density in heavy nuclei would 
rise cx A^/^^ but the nonlinear term rises faster, cx A^/'^, 
and slows down the A-dependence. Notice that the GLR- 
MQ equation with a quadratic nonlinear term describes 
only the onset of saturation. While Gn{x,Q'^) increases, 
the higher order terms become important. Eventually, at 
large A the gluon density is expected to reach a saturated 
value Gj^{x,Q^,b), which is independent of b (except at 
the nuclear periphery) . Thus, in the saturation limit the 
nuclear ratio of gluon densities depends on A as. 



K{x,Q' 



JdHG%{x,Q^,b) ^ ^_i/3 



AGn{x,Q^ 



(5) 
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The same phenomenon of saturation can be under- 
stood in the rest frame of the nucleus, which gives a more 
familiar interpretation having direct ties to the usual in- 
terpretation of shadowing in optics. In this frame gluon 
shadowing can be treated as Glauber shadowing for a 
glue-glue dipole [3] , or as a contribution of higher Fock 
components of the projectile containing gluons [s']. In 
this frame saturation corresponds to the black disk limit 
for the dipole-nucleus amplitude. 



III. MODEL INDEPENDENT UNITARITY 
BOUND 

There must be a bottom bound for the saturated den- 
sity, i.e. for the strength of gluon shadowing. Indeed, 
how many gluons are left in a heavy nucleus at small x 
after many fusions? Can it be as little as in one nucleon? 
As far as the nucleon is much smaller than the nucleus, 
this is obviously impossible. Two nucleons separated in 
impact parameter plane by a long distance ~ i?^ do not 
" talk to each other" even at small x [l6| . There should be 
some minimal amount of gluons at small x, correspond- 
ing to a maximal fusion rate at a given impact parameter. 
Assume that we have a one-dimensional row of nucle- 
ons. Intuitively, even at very low x I/ttinRa, where 
all gluons originated from different nucleons, overlap and 
fuse, the resulting gluon density cannot be smaller than 
in a single nucleon. Indeed, upon approaching the gluon 
density in a single nucleon, the further evolution in x de- 
scribed by Eq. [2] proceeds exactly like in a nucleon. So it 
is clear that the minimal amount of gluons in the nucleus 
should be proportional to its area, i.e. to A^^^, but the 
coefficient is important. 

Gluons originating from different nucleons and overlap- 
ping in longitudinal direction can fuse if their transverse 
separation is within the gluon interaction radius rg{x). 
The latter can be interpreted as the radius of the triple- 
Pomeron vertex. If the radius of gluon distribution in 
the nucleon is ris[{x), the gluon clouds of two nucleons 
overlap, if the distance between the nucleon centers does 

not exceed rjvAf = ^ r'jf{x) + rg{x). The mean number 

of nucleons whose gluons are able to fuse at impact pa- 
rameter b is. 



' NN 



TAib). 



(6) 



Assuming uncorrelated Poisson distribution, the prob- 
ability to find n nucleons in the tube is W„ = 
exp(— (n)) (n)"/n!. Maximal gluon shadowing corre- 
sponds to full fusion of gluons belonging to different nu- 
cleons in the tube, i.e. the collective gluon density of 
n > 1 overlapping nucleons equals to a single nucleon 
gluon density, nGNix,Q^) => Gn{x,Q'^)- The probabil- 
ity to have one or more nucleons in the tube, i.e. to have 
a saturated gluon density, is 1 — exp(— (n)). Correspond- 
ingly, the shadowed gluon density corrected for finite A 



has the form, 



GNix,Q^)^GN{x,Q^) { 1 



Hb)) 

A 



, (7) 



and the nuclear gluon density per unit of transverse area 
is 



GAix,Q^b) 



G{x,Q^ 



(8) 



' NN 



Correspondingly, the bottom bound for the nucleus-to- 
nucleon ratio of gluon densities reads. 



/ d%GA{x,Q\b) 
AGjv(x,Q2) 
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(9) 
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eff _ 



Jd'b{l - [1- 

cr'^J'^TAib)]^}. Here we maximized the denominator, 



where a'^J'^ — 
a"n ^ 

assuming fully saturated gluon density in a nucleon. For 
heavy nuclei this bound Eq. ^ can be approximated as 
Rj{x,Q^)>{R\/ATl^)^A-^/\ 

Now we should settle the values of Tg and tat which are 
functions of x and . Here we should rely on available 
experimental data. Gluons originated from two nucleons 
with different impact parameters overlap, if the mean 
distance squared, r'j^{x,Q'^), equals to the doubled mean 
radius squared of the gluon distribution in each proton. 
The latter was measured at HERA in reaction of virtual 
electroproduction of vector mesons. 



\rl{x,Q^) 



2Bei{x,Q^) 

2Bo{Q^)+4a'jp{Q^)ln{l/x). (10) 



The a:-independent term i?o ~ 4 GeV~^ is related to the 
proton formfactor; the effective slope of the Pomcron tra- 
jectory a'jp{Q^) varies with Q^. It is known to be as 
large as 0.25 GeV~^ in soft hadronic processes, and is es- 
sentially a result of unitarity saturation at small impact 
parameters [1, [l3| . Data on exclusive electroproduction 
of vector mesons show that at high the effective slope 
is much smaller, aip(Q2) « O.lGeV-2 (see more 

references in [3| ) ■ To be on the safe side (regarding model 
independence of the bottom bound) we use the maximal 



value 



= 0.25 GeV" 



The radius of the triple-Pomeron vertex, rg{x), is well 
measured in soft diffractive dissociation of protons to 
large invariant masses pp — > pX available with high 
statistics. The triple-Regge analysis of data 20] shows 
that this radius is compatible with zero, in any case is 
much smaller than r^- This observation gets a natu- 
ral explanation Q if the mean transverse momentum of 
gluons in the proton is large and controlled by a satura- 
tion scale Qsix) which increases as a function of energy 
[l3,[3, HH. Therefore, the radius rg{x) can only shrink, 
so we can safely neglect it compared to the nucleon radius 

TN- 
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With these values of the parameters the bottom uni- 
tarity bound, Eq. for gluon suppression in gold at 
= 1.69 GeV'^ is depicted in Fig. [1] The nuclear ratio 
calculated with EPS08 is also plotted in this figure 
and is fully in the shaded area below the unitarity bound. 

One can also trace the A-dependence of gluon shad- 
owing. In Fig. O we plotted the effective number of nu- 
cleons, Aeff = ARf{x,Q'^), suggested by the EPS08 
analysis at = 1.69 GeV^ and x = IQ-'^, 10"^. A^ff 



- Q^=1.69GeV^ 




A 

FIG. 2: Unitarity bounds for the effective number of nucleons, 
Aeff = ARf, as function of A calculated with Eq. ((Ojl at 
= 1.69 GeV^ for x = 10^" (top thin line) and x = 10"'' 
(bottom thin line). Solid thick curves are calculated with the 
EPS08 code Q. 

is quite small, especially at x = 10^® where for gold, as 
we already mentioned, Aef / = 1 • The A dependence also 
looks odd: A^ff never falls with A, it can only rise, as 
j^2/3 pj. fasi^gp^ Apparently, the predicted A-dependence 
at X = 10^® is unacceptable. Although it rises with A 
aX X = 10^**, but too slow, breaking the miitarity restric- 
tion. 

We also calculated the A-dependent unitarity bounds, 
Eq. ([5]), and plotted the results in Fig. Oat x — 10~^ and 
X — 10"^ by thin solid lines, top and bottom respectively. 
Again, the results of EPS08 shown by thick solid curves 
are quite below the bounds. 



at very small x (due to high-mass gluonic fluctuations), 
and is certainly not the case, say, at a: ~ I0~^ where the 
coherence length for gluons is shorter than the nuclear 
radius [225 . this range of x the bottom bound should 
be considerably lifted up. Actually, the onset of gluon 
shadowing starts only at a; < 10~^ Q. We exaggerated 
the strength of nuclear shadowing for the unitarity bound 
in order to keep it model independent. 

We also maximized the strength of gluon fusion, as- 
suming that the appearance of an extra gluon results 
in an immediate fusion, so that the gluon density re- 
mains unchanged. This corresponds to a black disk limit 
for pp elastic partial amplitude, or full saturation, while 
data show |23j that this regime is relevant only to central 
(6 « 0) pp collisions. 

Then, we maximized the glue-glue correlation radius, 
which we evaluated at rjv ~ 1 fm, while calculations 
on the lattice lead to a much smaller value of about 
0.2 — 0.3 fm [131 ■ There are also numerous evidences in 
data 0, |T3| indicating that gluons are located within 
small spots of radius ro ~ 0.3 fm around the valence 
quarks. Although, following the valence quarks, they 
are distributed in impact parameters with large radius 
Eq. (jlOp . the probability to overlap in impact parame- 
ters is significantly reduced. Indeed, for a given gluonic 
spot the mean number of others to overlap with is, 

{ng)^^rmTA), (11) 

where the mean nuclear thickness, {Ta) = \ ^ (PhT\{h) 
is about 1.5 fm^^ for heavy nuclei. Thus, transverse 
overlap of gluonic spots is quite small even for heavy 
nuclei, {ug) « 1.3. To improve the unitarity bound one 
can rely on Eq. ([5]), but with a different value of the 
effective cross section, cr^^ a^j^ = Ztttq. The new 
bound for gold at ^ 0.169 GeV^ Rf > 0.6, is about 
three times higher than the model independent one, 
depicted in Fig. [TJ 



V. GLUON SHADOWING FROM THE 
BALITSKY-KOVCHEGOV (BK) EQUATION. 



IV. MORE STRINGENT CONSTRAINTS 

Trying to be model independent we apparently exag- 
gerated some estimates. We assumed in ([7]), (|9]) that 
gluons in different nucleons completely overlap in longi- 
tudinal direction, i.e. that the coherence length is much 
longer than the nuclear size. This is questionable even 



It is interesting to compare the unitarity bound for 
gluon shadowing with expectations coming from contem- 
porary models of gluon saturation in nuclei. Although 
the GLR-MQ equation ([2]) describes qualitatively the 
phenomenon of gluon saturation and allows to take ex- 
plicitly into account the nuclear target, the correct equa- 
tion that governs the shadowing effects (BK equation [2^ 
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) has a more complicated form, namely, 
dN{r,Y;b) _ Cpas 



dY 



(12) 



2N[r',Y]h~ ]^{r-r')] - N (r,Y;b 



N {r,Y;b-^{f~f")j N (t^ - f' ,Y;b - 



where Y = ln(l/a;). This equation has a form differ- 
ent from the GLR-MQ, since the nuclear target enters 
only in the initial condition at some rapidity Y = Yq. 
This happens because the equation describes a simple 
process: a small dipole of size r decays into two dipoles 
of sizes r' and \r — f'\ with probability K(r,r') = 

Cp as 



y ^,2 , which does not depend on the tar- 
get. These two dipoles can interact either independently 
( the linear term of the equation), or simultaneously 
( the non-linear term of the equation). In the case 
of the nuclear target one can simplify Eq. taking 



Ra > r'andl 



So 6 becomes a parameter 



in the equation, and the ^-dependence is entirely deter- 
mined by the initial condition. In other words, there is no 
correlations between partons with different b. It should 
be noticed that an analysis of HERA data based on the 
BK equation was performed in the entire kinematic range 
of Q2 from Q2 = q.I GeV^ to = 100 GeV^ with ex- 
cellent x^/d-o.f < 1 

The key property of the BK equation (as well as of 
the GLR-MQ one) is the a ppe arence of the new dimen- 
sional scale Qs{x; b) [Tsl. [2ll. l27j . For further discussion it 
is worth mentioning that the saturation momentum for 
nuclei equals to 



It turns out that inside the saturation domain and even 
outside, the solution of the BK equation N{r, Y; b) obeys 
geometrical scaling |28|], i.e. depends only on one vari- 
able, rQs{x,b). In Ref. [131 it was found that in the 
saturation domain it has the form. 



N{r, Y;b) = 1-C cxp 



1 



7c 



2x(7c 



(14) 



where z = In (r^ Ql{x; 6)/4), and xil) is the BFKL ker- 
nel (function K{r,r') in 7, i.e. anomalous dimension) 
representation; ^cr is the critical value of anomalous di- 
mension corresponding to the saturation momentum (see 
details e.g. in Ref. [13]). Eg. p4)) shows how approaching 
the unitarity boundary depends on the target, namely, 
through its saturation momentum. Unfortunately , the 
value of the coefhcient C in Eq. can be taken only 
from the numerical solution of the BK equation, and it 
was found in Ref. [131 to be equal to C = 1/e. Further, 



the ratio i?^ can be written in the form 



Rf{x,Q) - 



fd^r\^g{Q;r)\^NA{r,x;b) 
A Jd^r\'i'g{Q;r)\^Np{r;x;b) 



(15) 
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Q =1 .69 GeV' 
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FIG. 3: Gold to proton ratio of gluon densities at = 
1.69 GeV^ as function of saturation momentum squared Q^{x) 
(see Egs. I13land ll6p for two models for the dipole amplitude 
N{r,x;b): solid curve for the solution of the BK equation in 
the saturation domain (see Eq. (|14p ); and dashed curve for 
the phenomenological GBW model (see Eq. (HHJ). 



where |^'g(Q;r)p is the probabity to find a gg dipole in 
the probe. This wave function has been discussed in [^, 
HH. It turnes out that |*g(Q;r)|2 oc for r > 1/Q. 

For N we use Eq. (jl4p with the saturation momentum 
for a nucleus given by Eq. (jl3p , and for a nucleon we take 



{A;x-M = dH'Qs{N-x-b')TA{b) (13) with 



{N;x-b) = Qs{x)S{b) 



Sib) - (^] kJ'^'' 



(16) 



R 



R 



One can see that S{0) —1 and with R = = 
0.89 fm it gives correct electromagnetic radius of the 
proton. 

For N{r,Y;b) in Eq. ^ we take Eq. ^ with 
r^Qs{x,b)/4: > 1. Therefore, in Eq. (fTS]) we should inte- 
grate over r > Tmin where 



Q Qs {N- X- b) J ' 



maximum 



Qs {A]X]b) 



(17) 



The predicted ratio is shown in Fig. [3] as a function 
of Qs = Qs{x) in Eq. ^ and Eq. while in Fig. H 

we show the prediction for the saturation momentum as 
function of ln(l/a;) as it was calculated in [30i|). 

Notice that our calculations are valid only for > 
1/Qg, a condition which is broken at small Qs. This is 
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5 7.5 10 12.5 15 
ln(1/x) 

FIG. 4: Saturation momentum Qs(x) for the proton versus 
ln(l/a;) . Curves 1 and 2 describe Qs with LO and NLO 
BFKL kernels (see Ref. [sj) respectively. Curve 3 presents 
the saturation momentum that stems from the modified BK 
equation of 30]. Curve 4 is the GBW phenomenological sat- 
uration momentum. 

the reason for the odd effect of antishadowing which one 
can see in Fig. [3] at small Qs- 

We also plot in Fig. [3] the prediction for from 
the model for the dipole amplitude N{r,x;b)^ inspired 
by the Golec-Biernat and Wusthoff (GBW) parametriza- 
tion, namely, 

Nir,Y;b) = 1 - exp (-r^ g2(a;; 6)/4) 

This model describes quite well the DIS data of HERA 
and comparison with this model can illustrate how close 
are our theoretical estimates to reality. 

VI. SUMMARY 

• Although at small x gluons originated from differ- 
ent nucleons in a nucleus overlap in longitudinal 
direction, most of them are still well separated in 



the transverse plane, therefore cannot fuse. For 
this reason the gluon density in nuclei cannot drop 
at small x below a certain bottom bound which 
we evaluate assuming maximal strength of gluon 
fusion. We found a model-independent unitarity 
bound for gluon shadowing. Examples for the nu- 
clear ratio of gold-to-proton and the A-dependence 
of the bound are shown in Figs. [T] and [5] respec- 
tively. 

• The recent results of the global fit of nuclear PDFs 
[l| lead to a gluon suppression which is orders 
of magnitude below the unitarity bound. More- 
over, the amount of effective nucleons predicted by 
EPS08 has an odd A-dependence (rising too slow 
or even falling with A) , and is also much below the 
unitarity limit, as is shown in Fig. [21 

The source of the problem is rather obvious. Try- 
ing to access minimal Bjorken x for nuclear par- 
tons, the EPS08 analysis included RHIC data on 
high pt hadron production at forward rapidities. 
However, these data were analyzed in [l[ assum- 
ing that the mechanism of suppression is related 
only to a modification of the nuclear PDFs, while 
at least two different mechanisms contribute to the 
observed nuclear effects. 

• Since maximal gluon shadowing is expected in the 
regime of gluon saturation, we performed model 
calculations based on the Balitsky-Kovchegov equa- 
tion and found that the results presented in Fig. [3] 
are well above the unitarity bound. 
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